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r) Abstract. We analyse the multifield behaviour in D-brane inflation when contributions from 



^ 



o 



the bulk are taken into account. For this purpose, we study a large number of realisations of 



Qj the potential; we find the nature of the inflationary trajectory to be very consistent despite the 



complex construction. Inflation is always canonical and occurs in the vicinity of an inflection 



?-H point. Extending the transport method to non-slow-roll and to calculate the running, we 

c/^ obtain distributions for observables. The spectral index is typically blue and the running 

i__i positive, putting the model under moderate pressure from WMAP7 constraints. The local 

/nl and tensor-to-scalar ratio are typically unobservably small, though we find approximately 

^ 0.5% of realisations to give observably large local /nl- Approximating the potential as sum- 

CN separable, we are able to give fully analytic explanations for the trends in observed behaviour. 

^j Finally we find the model suffers from the persistence of isocurvature perturbations, which 

(^ can be expected to cause further evolution of adiabatic perturbations after infiation. We 

*y^ argue this is a typical problem for models of multifield inflation involving inflection points 

<^ and renders models of this type technically unpredictive without a description of reheating. 
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1 Motivation 

Inflation (for recent reviews see Ref. [1]) is widely viewed as the most elegant paradigm to 
understand the very early universe, but despite being a simple set-up, it is hard to fully 
describe it within an ultra-violet complete theory. The main reason for this is that its 
dynamics are highly sensitive to Planck-scale physics. Generally speaking, Planck-suppressed 
contributions arise by integrating out heavy fields, which are present as extra degrees of 
freedom necessary for the ultra-violet completion of a theory. Unless protected by some 
specific symmetry, the inflaton will couple to these heavy fields, which results in a radiative 
instability of its mass and of the fiatness of the potential. It is therefore of enormous interest 
to try to study inflation in an effective action that takes into account contributions of high 
energy physics. 



In string theory, the important degrees of freedom to take into consideration are the 
heavy moduh that arise from stabiUzed compactifications of the extra dimensions. To com- 
pute the detailed contributions that their couphng to the inflaton induces in the effective 
action requires full knowledge of the stabilized compactification, which is rarely possible. In 
this sense, it is important to identify and work with string set-ups that present a sufficient 
level of computability. This is the main motivation to look at inflation arising from the 
dynamics of D3-branes in warped throats. This scenario is not expected to be generic but 
it allows some major simplifications that make the task of building the effective action more 
achievable. 

The flux compactification causes warping of the manifold, giving rise to regions of the 
bulk with warped throats. Inflation can occur in this scenario when a D3-brane, correspond- 
ing to our four-dimensional space-time, is Coulomb attracted to an anti-D3-brane that sits 
at the infra-red tip of the throat, where it minimizes its energy. Inflation is driven by the 
dynamics of the D3-brane, and it behaves like a multi scalar field system, where the fields 
can be viewed as the physical coordinates separating the branes. The infiationary epoch ends 
when the branes collide and annihilate. 

What is special about these warped throat regions is that they can be approximated 
by a finite region of a non-compact conifold geometry, for which the metric and background 
fluxes are well known. This finite segment is then glued to the compact bulk at some ultra- 
violet scale. Corrections to this non-compact approximation will arise from the effect of fields 
on the throat, like stabilized moduli of the compact bulk, and are clearly examples of the 
Planck-suppressed contributions to the effective action mentioned above. 

Even working in this 'simple' set-up, the effective action cannot be fully computed. At 
most, it is possible to calculate the form of the Planck-suppressed contributions, but they 
will come necessarily with an unknown Wilson coefficient. For this reason, to study the 
inflationary potential that arises from the D-brane scenario one needs to sample over a lot of 
realisations of random Wilson coefficients. Observational predictions need to be understood 
in the light of this statistical nature.^ 

Another issue to keep in mind is that the conifold approximation for the throat does not 
hold at the infra-red tip. This implies that, if using such an approximation, the dynamics can 
only be analysed before the tip, ignoring everything that occurs around and just before the 
collision of the branes. One can hope that this regime will not significantly affect the curvature 
perturbation (", otherwise the predictions made are useless for comparisons with observations 
at decoupling time or after. To check if this is a reasonable assumption, it is necessary to 
keep track of isocurvature modes as they transfer power to the curvature perturbations. If 
they have not completely decayed by the end of the analysis, the curvature perturbation will 
continue to evolve as the brane moves into the tip. The inflationary trajectory is said to 
have not yet reached its 'adiabatic limit'. The observable properties of C, will then depend 
on unknown details of the tip, including reheating, making the model incomplete. 

In this paper, our aim is to exhaustively study the possible D-brane dynamics above 
the tip and understand the consequent inflationary behaviour. For this purpose, we follow 
the most sophisticated set-up in the literature [2]. This includes the sampling over Wilson 



^This issue is analogous to the challenge of making predictions in models of the string landscape. In order 
to compare predictions from this model with observation we must assume we are a typical observer. Of what 
class of observer we are typical is however a very difficult question to address and brings with it an inherent 
measure problem. We make no attempt to address this interesting challenge here, but it should be noted that 
a resolution of this problem will add a weighting to the distributions we present in this work. 



coefficients and assumptions like the ones described above. The throat is approximated by 
a conifold parameterized by one radial and five angular directions that can be effectively 
viewed as the scalar fields driving inflation; this approximation only holds for a fixed range 
in the radial direction and it can be shown that within it DBI effects can be neglected. The 
framework and details of the construction of the potential are reviewed in §2. 

Since the potential that describes the motion of the D-brane is sensitive to all 6 co- 
ordinates, multifield effects have a profound impact on the dynamics of the inflationary 
trajectory and consequent curvature perturbations ("• To compute observables within this 
multifield super horizon dynamics, we use the transport equations method originally intro- 
duced in Refs. [3, 4] along with an extension to non-slow-roll that will be described fully in a 
forthcoming paper [5] . This technique is a realisation of the separate universe assumption in 
which the values of correlation functions of C and their tilts can be directly transported from 
horizon-crossing to the desired time of evaluation, as described in §4. Since this is framed 
in terms of ordinary differential equations, it allows for a clean and efficient numerical im- 
plementation as well as a means of explicitly keeping track of the evolution of isocurvature 
modes. 

As a consequence of our analysis, we obtain the probability of getting inflation, found 
to be in agreement with Ref. [2], and the statistical distributions of observables predicted 
by our set-up. But the most interesting results of this work come from our understanding 
of the peculiarities of the multifield dynamics of the D-brane potential. By looking in detail 
at the resulting inflationary trajectories, we are able to map trends in the distributions of 
observables to generic features of the potential. Moreover, we can see if these features allow 
the trajectories to reach their adiabatic limit. The outcome of our work is then more than 
the computation of predictions, being an investigation on the limitations of our set-up as a 
predictive and useful toy model for the D-brane scenario. 

2 D-brane inflation 

2.1 D-branes in a warped throat 

As mentioned earlier, in the D-brane inflation scenario, our Universe lives in a D3-brane and 
experiences inflation due to its dynamics in a warped throat region of a stabilized compact 
space; the D3-brane is attracted by an anti-D3-brane that sits at the tip of the throat where 
its energy is minimized. When the two branes collide, they annihilate, inflation ends and 
reheating occurs [6, 7]. 

The throat region, excluding the tip, can be well approximated by a non-compact coni- 
fold geometry; the conical singularity that would arise at the tip is smoothed by fluxes such 
that the radial coordinate at this point is finite. For the purpose of this work, only the region 
above the tip will be considered. In this case, and ignoring logarithmic corrections to the 
warp factor, the background geometry can be approximated by: 

ds' '^ 



") gijd^'d(t^ = (^^ {dr^ + r^dslui), (2.1) 



where r is the radial conical coordinate and T^'^ is the coset space {SU{2) x SU{2)) /U{1) 
that describes the angular directions of the cone. The radius R is approximately ruv; the 
coordinate at which the throat is glued to the compact bulk, as illustrated in Fig. 1. In 
agreement with Ref. [2] we use the value ruv = 1 throughout this work. It is useful to define 
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Figure 1. The non-compact conifold approximation for the warped throat. This approximation 
holds between x = 1, where the throat is glued to the compact bulk, and x = 0.02, where the tip that 
cannot be described by our approximated geometry starts. 

a rescaled radial coordinate as x = r/ruv that in the cone region is always ^ x < 1. 
To ensure that the non-compact approximation always holds, we restrict our analysis to the 
regime comfortably above the tip, where 0.02 < x < 1. The value x = 0.02 was chosen in 
agreement with Ref. [2]. 

The T^'^ space is parameterized by 5 angles ^ = {9i,92,^i,^2,i^} where < 6'i < vr, 
<92<TT,0 < ifi <2Tr, <ip2 <2tt and < -0 < 47r. 



Throughout this paper we will use units M, 



Ref. [2] , we use the value ruv 
Lagrangian like 

C = a' 



py — SttG = 1. Also, in agreement with 
1. In this throat the D3-brane experiences a DBI inflationary 
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where a is the scale factor, T3 is a constant representing the brane tension and, within 
the approximation where the logarithmic corrections to the warp factor can be ignored, 
T((/)) = Tsx^. The value of the warp factor at the tip is determined by the parameter ao 
such that T((/))|tip = TsOo- Some physical arguments concerning the consistency of the set-up 
enforce a limit on how small T{(j)) can get [8]; in this work, following Ref. [2], we use the 
values T3 = 10-2 and oq = lO'^. 

As mentioned in Ref. [2], for our specific realisations of the D-brane action, the brane 
velocity is always very small compared to T{(f)), making 
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This is equivalent to saying that DBI effects are negligible, as we can rewrite the Lagrangian 
as 



C 



a 



Tsgij^'^ - V{(l)) 



(2.4) 



and identify the canonical kinetic term rescaled by the constant T3. 

The fact that this simplification can be made is related not only to the choice of T3 
and ao but also to the fact that our analysis only includes the throat region above the tip. 



Condition (2.3) breaks if T((p) gets very small and, in fact, T((j)) decreases with x. Possible 
DBI inflation in this regime can have strong repercussions in the value of observables at the 
end of inflation; the inclusion of the whole throat in the computation of perturbations is then 
potentially very interesting but beyond the scope of this work. 

2.2 D-brane potential 

In the simplest form of this scenario, the potential that induces inflation has two contri- 
butions [7]. First, there is the Coulomb interaction between the pair of branes, which is a 
multipole expansion where high-multipole terms are suppressed by powers of oq . The leading 

terms are: 
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where the parameter Dq = 2T30^ determines the overall scale of inflation. 

Second, the coupling to the curvature induces, at leading order, a mass term like: 

Vm = ^/x^ (2.6) 

where the scale /i^ = DoT^r^jy/Mpy 

In this basic picture, inflation would actually be single-field, with the inflaton being 
the radial separation between brane and anti-brane, x. However, if the dynamics was deter- 
mined by these two terms only, sufficient inflation couldn't be achieved [6]. Such a potential 
has a single inflection point; this feature will be shown in the next sections to have strong 
consequences for the phenomenology of the full model. 

This simplified picture ignores some important contributions to the potential experi- 
enced by the D-brane. One needs to take into account that the throat is finite and glued 
to a compact manifold, and as such, moduli stabilization from the bulk will necessarily have 
an impact on the throat geometry. These contributions can be viewed as corrections to the 
non-compact approximation, and will be denoted by Vbuik- 

Ideally, one would like to have the full knowledge of the 4-dimensional potential induced 
by the compactification flux on the brane dynamics, but this is not possible to achieve for a 
general Calabi-Yau bulk. However, it is known [8] that such a potential, in the non-compact 
background of the conifold, respects the Laplace equation: 

Wbuik = 0. (2.7) 

Since we know completely the geometry of the conifold, this equation can be explicitly solved. 
We refer to these contributions, following the notation of Ref. [8], as the homogeneous con- 
tributions to Fbulk- 

Deviations from this expression, which holds for the non-compact background, can be 
obtained by allowing a source from the bulk. In this case, the Poisson equation looks like [8] 

Wbuik = g|A|2, (2.8) 

where gs is the string coupling constant and A is proportional to the imaginary anti-self-dual 
three-form fiux from the bulk. We will refer to these contributions as the inhomogeneous 
contributions to Vjauik- To solve this equation, a simplification can be used. Since these 
contributions are perturbations to the non-compact approximation, they can be assumed. 



up to a good approximation, to have the same structure as the homogeneous contributions. 
So the idea is to express them as an expansion of harmonic terms from the homogeneous 
solution. In other words, the solutions to the Laplace equation dictate the structure of the 
bulk contribution to the potential. 

Homogeneous contributions: The Laplace equation (2.7) for our non-compact coni- 
fold can be written in the form of the expansion [9] : 



Horn bulk(2;, ^) = / ^ CLMX^^^'^YLMi^) (2. 



LM 

where Clm are constant coefficients, YlmI^) are the angular eigenfunctions of the Laplacian 
of the T^'^ space and the subscripts L = {li, I2, R} and M = {mi, 7712} represent the quantum 
numbers under the T^'^ isometries. The powers A(L) are related to the eigenvalues of the 
Laplacian and are given by 

A(L) = -2 + ^6/i(/i + 1) + 6/2(^2 + 1) - 3R/4 + 4. (2.10) 

The magnitudes of Clm are highly dependent on details of specific compactifications, 
so they need to be considered unknown parameters. Using the scale //^, considerations of 
Ref. [10] suggest that Clm ~ C'(l)i so a way to deal with this lack of knowledge is to scan 
randomly over values in this range. To take the leading contributions of this term, one needs 
to consider the lower values of A(L). The maximum value desired for A(L) determines the 
truncation of the summation. 

Inhomogeneous CONTRIBUTIONS: To solve Eq. (2.8) as an expansion of the type of 
Eq. (2.9), one needs to identify the radial scaling of the flux A in terms of the quantum 
numbers L and M of T^'^. It is possible to classify the flux in 3 different series, I, II and III, 
regarding their different radial scaling [8]. The radial scaling of |Ap is given by ^ 

A(Lq,, -L/3)inhom bulk = Aq(Lq,) + /^^{Lp) — 4 (2-11) 

where a and (3 run over the 3 different series I, II and III and 



/\i{L) = -1 + V6/i(/i + 1) + 6/2(/2 + 1) - 3i?/4 + 4, (2.12) 



A//(L) = ^/Ui{li + 1) + 6^2(^2 + 1) - 3i?/4 + 4, (2.13) 

A//7(L) = 1 + ^6/i(/i + 1) + 6/2(^2 + 1) - 3ii/4 + 4. (2.14) 

It is then possible to write the inhomogeneous contributions as: 

l^nhom bulk(x, ^) = / Y. C'i,M.i^M,X^(^-^'3)yL,Mj*)yL,M,(*). (2.15) 

LaMcLpMfj 

To write this expression with the same structure as Eq. (2.9), the angular part needs to be 
expanded in terms of Ylm of T^'^ as 

YL^M^{^)YL,M,m = Y.A^pYlm{^) (2.16) 

LM 
^Some contractions of flux series vanish following the considerations of Ref. [8]. 




0.8 ---^0 
Figure 2. A typical realisation of the D-brane potential with only one angular direction, ■0, active. 



such that finally, 



Mnhom bulk(^, ^) = / Yl Y.^LMx'^^''-'''^^A^(sYLM{'i'). 



(2.17) 
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The constants Clm correspond to the random parameters associated to each Ylm from 
the homogeneous contribution, and, just as in that case, the maximum value desired for 
A(Lq, Lp) determines the truncation of the summation. 

The total potential experienced by the D-brane in the throat is then 



V{x, ^l) = Vc + Vm + Vhom bulk + tinhorn bulk 
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+ /i'5]Cz.MX^WyiM(M/) + / Yl Y^lmx'^'^'^-''''^A^pYlm{^) 
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A specific realisation of this potential, with only one angular direction being taken into 
account, is shown in Fig. 2. 



3 Experimental procedure 

Having established that the brane potential necessarily has a level of randomness to be able to 
encompass complex contributions, it is important to construct a useful sample of realisations 
for the study of the emergent inflationary behaviour. 

The first thing to specify is the maximum values of A(L) and A(Lq, L^) in the potential. 
Since the mass term has power ex x^, it makes sense to include at least all terms with A < 2. 
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Figure 3. Distributions for tiie scalar spectral index for the model with 2 active fields. The distri- 
butions were taken with different values of Q and show very similar behaviours (Amax = 3 in all 
distributions). 



For computational reasons, and in accordance with Ref. [2], we looked at potentials with 
A < 3 and with A < \/28— 3/2. This corresponds to a total of 121 and 390 independent terms 
in the potential, respectively.'^ The values taken by A are: 1, 3/2, 2, 5/2, v28 — 5/2, 3, v28 — 
2, 7/2, \/28 - 3/2. 

The second thing to decide is how to generate the random Clm coefficients. Following 
Ref. [2], we define Clm = QClm such that Clm is a distribution with unit variance, en- 
capsulating the information on the distribution, and Q is the root mean square size of Clm-, 
encapsulating the information on its magnitude. 

Drawing conclusion on the predictions of the model would be problematic if the infla- 
tionary behaviour emerging from our sample was dependent on the type of distribution of 
Clm- Fortunately, as shown in Ref. [2], this is not the case. In this work, we use a Gaussian 
distributed Clm- 

Regarding the choice of Q, a similar argument could be invoked. As mentioned in the 
previous section Q ~ 0(1), but the probability of inflation is very sensitive to its precise 
value [2]. If the emergent behaviour, given that inflation occurs, was also sensitive to this 
choice of Q, it would be hard to make general predictions. To ensure that no issue would 
arise from this effect, we tested the dependence of the inflationary phenomenology on the 
choice of Q and found it to be independent. This is demonstrated in Fig. 3, where, as an 
illustrative example, the distribution of the spectral index for a two-fleld sample is plotted 
for flve different values of Q. 

The last thing to flx before constructing the sample is the choice of initial conditions, for 
which we follow Ref. [2] precisely. Since the potential is statistically invariant under angular 
translation in T^'^, generating multiple realisations starting always at the same angular 
coordinate automatically encompasses the statistical effect of varying initial conditions within 



^Note that these numbers differ slightly from those of Ref. [2]. The origin of this discrepancy lies in details 
of the expansions performed on Vbuik- 



a single realization. We hence consider just one initial condition per realization, arbitrarily 
taken as the angular coordinates being ^o = {Ijljl)!)!}- Regarding the initial radial 
direction, and following arguments from Ref. [2], we chose xq = 0.9. We set all the initial 
velocities to zero, xq = ^o = 0, leaving the study of the possible impact of initial velocities 
on inflationary phenomenology for future work. 

We can now present our experimental procedure for the building of a statistical sample 
of inflationary trajectories. This procedure was first used in Ref. [11] and more recently in 
Refs. [2, 12, 13]: 



1. Generate a random potential V{x, ^) starting at (xq, ^o) 
to find the field trajectory. 



(0.9, 1, 1, 1, 1, 1) and evolve 



2. If the model gets stuck in eternal inflation, i.e. does not reach x = 0.02, reject. 

3. If the brane gets ejected from the throat, i.e. x gets larger than 1, reject. 

4. Once the brane has reached x = 0.02, if the number of e-folds of inflation N < 60, 
reject, as insufficient inflation occurred, otherwise calculate observables. 

5. Repeat steps 1-4 many times to obtain a statistical sample. 

Several sets of samples were generated by changing some parameter, either AmaXi Q or the 
number of active flelds. When mentioning two-field samples, we are referring to a model with 
the radial direction and one angular direction active. When we discuss more active fields, we 
are referring to models with additional angular directions. 

The probability of achieving successful infiation in this set-up is in agreement with 
Ref. [2]. If Q is too small, most of the trajectories do not produce 60 e-folds of infiation; 
if Q is too large, most of the trajectories lead to an ejection of the brane. The optimum 
value of Q lies between these two regimes. As the number of active fields increases, just as 
when Amax increases, P{N > 60) becomes more sensitive to the choice of Q, and decreases 
slightly. Illustrative values are shown in Table 1. It is interesting to note that even for the 
optimum Q, more than half of the rejections are due to the brane being ejected from the 
throat, in agreement with Ref. [2]. 
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Table 1. Probability of getting successful trajectories {N > 60) as a function of the number of active 
fields and Amax- 

Regarding the value of Amax, we noticed that the inflationary phenomenology does 
not depend strongly on the chosen value. An illustrative example is shown in Fig. 4, for 
the spectral index of models with two and three flelds. For this reason, in what follows, we 
concentrate on distributions with Amax = 3. 

In the next section, we present the techniques used for the computation of observables. 



4 Computing the curvature perturbation 

A central feature of any model of inflation with more than one active scalar fleld is that 
the primordial curvature perturbation (^ evolves on superhorizon scales. To compute this. 
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Figure 4. Distributions for the scalar spectral index for the model with 2 active fields, left, and 3 
active fields, right. The distributions were taken with Amax = 3 and Amax — v 28 — 3/2 « 3.8 and 
show very similar behaviours ((5=1.4 in all distributions). 

almost all methods in the literature to date make use of some variant of the separate universe 
assumption [14-16]. The idea is to understand the final curvature perturbation as the result 
of the scatter of a collection of equal size space-time patches, within which the field values 
at a given time are constant. Let us make this more precise. 

4.1 The separate universe assumption is a geometrical optics approximation 

The separate universe assumption states that, when smoothed on some physical scale L 
much larger than the horizon scale, the average evolution of each L-sized patch can be com- 
puted using the background equations of motion and initial conditions taken from smoothed 
quantities local to the patch. The evolution of (" can be understood as the variation in the 
expansion of these patches. Since each smoothed patch corresponds to a position in phase 
space, the evolution of (^ can be determined by the evolution of an ensemble of points in the 
classical phase space. These points are subject to the laws of statistical physics and hence 
evolve according to the Liouville equation. 

Under the separate universe assumption, interactions between patches are negligible and 
therefore all that is required is a mapping of the initial conditions to a final state. The final 
distribution in phase space can then be viewed as an image of the initial conditions. This 
mapping simply follows a flow generated by the background theory and can be calculated in 
precisely the same way as geometrical optics enables us to calculate the image generated by 
a source of light rays. This optical description was made precise in Ref. [17], which we briefly 
summarise in this section, referring the reader to Ref. [17] for a more detailed discussion. 

4.2 Transport equations 

Since slow-roll approximations are not valid in general, we are required to work in a 2A^F-phase 
space. This consists of A'^p flelds (j)i as well as their momenta pi = (p^, where primes represent 
differentiation with respect to the number of e-folds N. The flelds <j)i and pi are treated on 
an equal footing, so from now on we will denote a point in phase space by ipa = {</'«, Pi} 
where a runs from 1 to 2A'^p. 

In canonical models of inflation, if we set initial conditions near horizon-crossing, the 
initial distribution of fleld perturbations will be close to Gaussian [18]. Furthermore, typical 
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spacing between arbitrarily selected members of the ensemble is of order the quantum scatter. 
It follows that the trajectories traversed by the ensemble trace out a narrowly-collimated 
spray or 'bundle' of rays in phase space with an initial Gaussian distribution. This scenario 
is well studied in the optics literature since many lasers have this characteristic. 

Cross-sections within the bundle of trajectories may be focused, sheared or rotated by 
refraction. It is ultimately through these distortions that any evolution in (^ occurs. To 
describe these distortions quantitatively, it is only necessary to know how some basis which 
spans the cross-section is transported from slice to slice. Denoting the difference between 
two field values at equal-time positions x and x + r by (5c/9a(r), we have an appropriate basis. 
This basis evolves along the beam as [17] 

^'^ = Ual3M^)]S^I3{r) + -Uaf3'y[v{^)]5ip(3{r)5ip-y{r) + ■■■ . (4.1) 

where Uq,^ = dpUa is the expansion tensor defined as the derivative with respect to the fields 
of the background fiow iiQ,[99(x)] = v'a(x) and similarly Ua/s-y = d^Ua/s- For clarity, we will 
drop the explicit <f{x) dependence from now on. 

The expansion tensor can be decomposed as a dilation 9 = tr Uap , a traceless symmetric 
shear (Tq^, and an antisymmetric twist coap, 

a 

Ual3 = ITf^a/S + CTa/S + ^a/S, (4.2) 

Dilation describes a rigid, isotropic rescaling of 5(po, by 1 + 0, representing a global tendency 
of the light rays to focus or defocus. The shear ctjj represents a tendency for some light rays 
within the beam to propagate faster than others. The twist LOa/3 describes a tendency of 
neighbouring trajectories to braid around each other. 

The observables of interest, like the power spectrum of C, its spectral index Ug, the local 
non-Gaussianity parameter /nl, etc, are related to the correlators (CC) and (CCO- So to 
compute the evolution of these quantities we need to know the evolution of the correlators 
of 5(pa- The full set of basis vectors contains all information required to determine the 
evolution of the bundle, encoded in Eq.(4.1) by the u-tensors. To obtain transport equations 
for the correlation functions simply requires reorganisation of this information. As was shown 
in Refs [3, 4, 17] this can be done in a number of ways. A particularly quick method is to 
acknowledge that provided the perturbations can be treated classically, we expect d{0) /dN = 
{dO/dN) for any quantity O. We can therefore immediately arrive at expressions for the 
two-point and three-point functions. Writing the two-point function as Ti^p = {dfa^ipjs), 
Eq. (4.1) implies 

"m' " \^^^^ ^ '^'^"^f^/ " """T^T^ + ^/37^7« + [> 3 p.f.] (4.3) 

Similarly, writing the three-point function as aa/37 = {SipaSipisSip-^), we get 



dN 



Uaxaxfi'y + UaXfi^xis^ti-y + cyclic (« -)> /3 -^ 7) + [> 4 p.f.]. (4.4) 



This forms a coupled set of ordinary differential equations which can in principle be 
extended to any n-point correlation function (see Ref. [19] for an implementation of this 
technique for the trispectrum). In the context of this work we only care about the power 
spectrum and bispectrum of C, so this set of equations encodes all the information we need 
for the understanding of superhorizon evolution of our observables. 
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4.3 Gauge tranformations 

Having seen how to compute the evolution of the field perturbations, now we need to relate 
these to the primordial curvature perturbation (^. This can be done by realizing that the 
curvature perturbation ( evaluated at some time t = tc is equivalent on large scales to the 
perturbation of the number of e-foldings N{tc, t*, x) from an initial flat hypersurface at t = t,,,, 
to a final uniform-density hypersurface at t = tc [20], 

C(tc, x) ~ 6N{t^, U,x) = N{tc, U,x)- N{tc, Q (4.5) 

where A^(tc,t*) = J^ Hdt. In the transport method, the hypersurfaces aX t = t^ and t = tc 
are chosen to be infinitesimally separated. Expanding 5N in terms of the initial field pertur- 
bations to second order, one obtains 

C(tc, x) = 6N{tc, t„x) = N,c, 8^1 + -N,a(3 {S^*Jip} - {5^*M})), (4.6) 

where repeated indices should be summed over, and N,a, N^^p represent first and second 
derivatives of the number of e- folds with respect to the fields v'^.^ The N derivatives are 
simply a gauge transformation from field perturbations to curvature perturbations. This 
gauge transformation only needs to be performed at the time of evaluation of Q. The fact 
that it does not need to be transported through superhorizon evolution is a great numerical 
advantage of this technique. 

It is straightforward to express the observables of interest at the time of evaluation within 
this formalism. The power spectrum is just related to the two-point correlation function of 
Q and can be obtained by [3, 4] 

P^c = A^,aA^,/3S„/3. (4.7) 

To obtain initial conditions at horizon-crossing, we set all the fields to be effectively massless. 
We tested the consistency of this assumption and found it to be valid in every realisation.^ 
The scalar spectral index, which expresses how the power spectrum changes with scale, is 

defined as 

d\iiPc_c_ 



dink 
where k^, is the pivot scale. This can be rewritten as 



(4.8) 



"^ Pec dink - NxN,j:,,' ^^-^^ 

where we have introduced the matrix n^^^ = dTi^tis/dlnk. Since the gauge-transformation 
factors N^a and the expansion tensors Uai3 are ^-independent (they depend only on the 
typical trajectory followed by the smoothed fields at each time) we can understand how 
the spectral index evolves on superhorizon scales [21]. The only necessary ingredient is a 
transport equation for the object UajS which is: 

dUaS d dTiaH 



^The subtraction of the correlation function in the second term is due to the fact that this covariance matrix 
corresponds to the contribution from disconnected diagrams which gives the vacuum energy. In Fourier space 
one only considers connected diagrams from the outset and thus the subtraction is already implicitly taken 
care of. 

^We would like to thank Sebastien Renaux-Petel for pointing out this possible issue. 
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In this work, for the first time, we apply an equivalent method to evaluate the running 
of the spectral index, which estimates how Us itself changes with scale. It is defined as 



d{ns 



1 



dink 



d 



k=k-t 



dink 



N,aN^pnaA _ N^aN^fsr^p 



P< 



a 



cc 



K - 1)' 



(4.11) 



where we have introduced the matrix r^^ = duap/dlnk. Again, to have the evolution of the 
running, we just need to specify the transport equation for Tq,^ which is: 



dvap d driap 



dN 



dink dN 



Uaxrxp + upxrxa- 



(4.12) 



The horizon-crossing initial condition for this expression can be obtained by calculating the 
derivative with respect to Ink of the horizon-crossing value of Uap given in Ref. [21].^ This 
derivative is to be evaluated at equal times. To do so, we first compute how the Uap associated 
to the pivot scale k^ changes with a variation in scale like A;* + 5lnk. This is 



dn, 



afS 



dink 



(-24c.+4eL.)5] 



aB\h.c. 



h.c. 



du, 



"7 V 1 ^'"^7 V 



dN 



^7/3 



dN 



-ya 



h.c. 



+ 2eh.c. {'Uay^'yl3+Ui3-y^^a)\h.c. 

(4.13) 
where e' = de/dN. 

The k* and k^, + 5lnk modes cross the horizon at different times but we are looking for 
the change at equal times. When compared at the same time the longer mode experiences 
slightly more evolution. For this reason, we then need to include an extra contribution that 
corresponds to this Ua/^ displacement. Remembering that at horizon crossing dink r^ dN, 
this is just duap/dN . The total expression is then 



'^a^lh.c. 



(-2eL.+4eL.)S 



Yi.cJ^aPWi.c. 






dN 



dN 



■ya 



(4.14) 



The local non-gaussianity parameter /nl is defined as [22] 

5 i?cCC 



h.c. 
c- 



/nl 



18 {Pcd' 



(4.15) 



where -B^^c; is the bispectrum of curvature perturbations, related to the three-point function 
of C- It is useful to decompose it as B(^(^(^ = 5(^f(i + ^CCC2) where 



and 



S, 



CCC2 



B(^C,C.l = N^a N,i3 N,^ aal3-y, 

3 



(4.16) 



(4.17) 



Eq. (4.16) is the intrinsic non-linearity among the fields, while Eq. (4.17) encodes the non- 
Gaussianity resulting from the gauge transformation to (" [3, 4]. 

Given these expressions for the observables, it is now only necessary to specify the gauge 
transformation expressed by the N derivatives. Following the procedure developed in Ref. [5] , 

®We would like to thank David Seery for discussions on this topic. 
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the expansion of Eq. (4.6), can be understood in two steps. First, 6N can be written as an 
expansion in terms of 5p, where 6 refers to a change from an initial flat hypersurface to a 
final uniform-density hypersurface:^ 

where p = 3H^. To obtain the derivatives of N as desired, one just needs to perturb each 
term of the above expansion in terms of the fields. The result is [5] 

^- = 2Hll3-er ^'-''^ 

^,0>0, - WTFTTT^ :v - ,rTA2^o :n2 : — + ^e , (4.21) 



N,P.P. = TTTT^T^ - .i.o ^. ( ^^^^ -66+ 12 ) , (4.22) 



^,0.. = - .^.:.r .. ( ^^^ - 2. + 6 ) , (4.23) 



2if2e(3-e) 4//4e2(3 _ g)2 \^ ^ 
2e(3-e) 4e2(3-e)2V e 

4F2e2(3-e)2V e 

where the derivatives are taken explicitly with respect to (pi and pi rather than the combined 
ipa, the objects V^i and V^ij refer to the derivatives of the potential with respect to the fields 

(pi and (pj and p[ = dpi/dN = d'^cpi/d'^N. 

4.4 The adiabatic limit 

This section has been dedicated to the understanding of superhorizon evolution of the cur- 
vature perturbations. What is yet to be addressed is at what point C, ceases to evolve. This 
happens when the trajectory has effectively reached its adiabatic limit, i.e. the model has 
become effectively single-field (see Ref. [23] for early discussions on this topic). 

One needs to ensure that the adiabatic limit is reached before the time of computation 
of observables, which generally is taken to be at the end of inflation. If isocurvature modes 
are still present at that point, the curvature perturbations will continue to evolve through 
an epoch of reheating, for which there is no precise knowledge. This would mean that the 
conclusions reached at the time of estimation would necessarily be incomplete, and the model 
not predictive. 

In the D-brane case, the situation is even more delicate, since the validity of the approx- 
imations of the set-up breaks down as one approaches the tip of the conifold, or, for practical 
purposes, at x = 0.02. If ( continues to evolve after that point, it will suffer changes due not 
only to reheating, but also an unknown background geometry and potential. 

Formally, the approach to an adiabatic limit can be understood via the parameter 
9 = ti Uq,j3 of Eq. (4.2) which describes the tendency for the bundle of trajectories to focus 
or dilate. The factor by which the bundle's cross-section has grown (or decayed) from its 
horizon-crossing value at Nq to the time N is given by: 



G(A^, iVo) = exp 



6{n)dn 

No 



(4.24) 



''We particularly thank David Mulryne for this result and discussions around this topic. 
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In a slow-roll two-field case, where the phase space is simply (pi, when — )• 0, or equivalently 
6 — )■ — oo, the trajectory is becoming effectively single-field [17]. However, since we are 
working with more than two fields and in a phase space composed by (pi and pi the situation 
is more subtle, and the relation not so direct [5, 17]. 

In the full phase space, it is instructive to understand how the bundles of (pi and pi 
individually behave. In particular, let us identify pi = Pi^{(pi) + Sj, where pf^ is the slow-roll 
attractor for the momenta pi and a function of (pi only and Sj are the momenta isocurvature 
modes. With these new variables, we can write [5] 

9 = e'^''{(P,) + 9'{si) (4.25) 

where 9 describes the dilation of the field bundle and 9^ describes the dilation of the s 
modes, i.e. how the momenta converge to their slow-roll attractor. We can also define Q^^ 
and 0* in analogy with Eq. (4.24). 

The interpretation of when the adiabatic limit has been adequately reached from the 
above quantities is not straightforward [5], but to infer when it has not been is often simple. 
For the purpose of this work it will be sufficient to know that if 0^^ > 1 then an adiabatic 
limit certainly has not been reached. 

5 Distributions for observables 

We now present the results obtained from computing the curvature perturbations for differ- 
ent samples of inflationary trajectories. The outcomes are distributions for the values of the 
cosmological parameters presented in the previous section: amplitude P^^, spectral index n^ 
and running of the spectral index of the scalar power spectrum, tensor-to-scalar ratio r, and 
local non-gaussianity parameter /nl- We compare these with constraints from observations; 
all constraint contours are 95% confidence limits using the WMAP 7 year data release com- 
bined with baryonic acoustic oscillations and supernova data [24-26]. We use as pivot the 
scale that crossed the horizon 55 e-folds before the end of infiation. 

5.1 Field number dependence 

Regarding the dependence of the infiationary behaviour on the number of active fields, we 
found that the effect of increasing the number of angular directions is negligible for the 
trajectories and consequent observables. Fig. 5 shows how the increase in angular directions 
changes the distributions of the tensor-to-scalar ratio and spectral index. We can identify 
a slight suppression in r as the number of active fields increases, as one would expect since 
more directions allow for more turns in field space which fuel the scalar power spectrum. The 
spectral index, which consistently peaks at 1, gets a shift towards red values as the number 
of fields increases. This is not the consequence of a change in the inflationary trajectories 
and it can be readily understood, as will be shown shortly. 

Interestingly, although an increase in the number of flelds allows in principle a wider 
range of trajectories, we did not encounter any difference in the distribution for the total 
number of e-folds when varying Np, as can be seen in Fig. 6. 

In the remainder of this paper, we will concentrate mainly on the results for two active 
flelds for which our sample is constituted of 1086 inflationary trajectories, using the sample 
for the full six-field case, composed of 93 trajectories, for comparison. 
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Figure 5. Distributions for the tensor-to-scalar ratio r, left, and scalar spectral index Us, right. 
The distributions were taken for different numbers of active fields. It is possible to identify a small 
suppression in r and a tendency towards redder values of rig (keeping the peak at n^ = 1) as Np 
increases. 
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Figure 6. Number of e-folds of inflation for different number of active fields. We can see that this 
distribution is unchanged by varying Np. 

5.2 Observables 

The amplitude of the power spectrum, as can be seen in Fig. 7, consistently undergoes a 
superhorizon evolution that we refer to as 'caterpillar' shaped, i.e. non-monotonic. This is 
interesting as this sort of evolution is very rare in purely random generated potentials [13] and 
so indicative of a dynamical trait common to all of our inflation realisations. The next section 
will explore its origin. As can be seen in Fig. 8, the histogram of P^^ has a smooth maximum 
at around 10^^, in agreement with observations (the WMAP value is ~ 2.5 x 10^^ [25]). 
This is not surprising as the overall magnitude of the potential is determined by the scale 
/i , which in turn is set by our choice of the throat length ruv- The important fact is that 
the distribution does not sharply peak at a precise value of P^^, indicating that there is no 
fine-tuning issue around this parameter. 

The spectral index shows a much clearer peak at n^ = 1, as seen in Fig. 8. Actually, 
the spread in Ug is not well approximated by a Gaussian; instead, two different populations 
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Figure 7. Superhorizon evolution of the amplitude of the scalar power spectrum P^^ for 50 trajectories 
with 2 active fields, coloured at random. The evolutions show a consistent non-monotonic growth 
that we refer as 'caterpillar' shaped. It is impractical to show the full evolution of all the trajectories 
due to the large range in scales even for this reduced sample; for this reason the plot has been cut-off. 

can be identified, one with n^ > 1 and one with n^ < 1. For the two-field ensemble, these 
correspond to ~ 84% and ~ 16% of the trajectories, respectively. As seen previously, an 
increase in the number of active fields enhances the number of red trajectories; the respective 
ratios, for the six-field sample are then ~ 50% each. In the next section we will address the 
dynamical characteristics of these two populations and the changes with the number of fields. 

A very interesting result comes from the computation of the running of the spectral 
index. Fig. 8 shows that the running tends to be positive and that it can take large values. 
This outcome, as will be discussed shortly, is extremely constraining. One could think that, 
since the running can take large values, the value of the spectral index changes a lot with 
the choice of pivot scale. Our approach is to think that sampling over a large number of 
inflationary potentials reproduces the effect of sampling over different choices of pivot scales, 
such that the final distributions for different pivot scales are actually identical. We tested 
this assumption and found it to be the case. 

The tensor-to-scalar ratio is always extremely small, as it is related to the slow-roll 
parameter e that remains <C 1 throughout the calculation. This can be clearly seen in Fig. 9. 
Furthermore, the usual single-field result for the brane case would be the relation 

r = 16^ (5.1) 

which corresponds to the green line. We can see how multifield effects change this result, by 
weakening this expression to an inequality. 

The local non-gaussianity parameter /nLj as can be seen in Fig. 8, is almost always too 
small to possibly be detected by any anticipated observation. In the full two-field ensemble 
only five trajectories yielded values of |/nl| > 1- Although these are highly unlikely, any case 
presenting interesting observational signatures can be informative in its own right; we leave a 
detailed analysis of these cases for future work. When plotted against Ug, this parameter also 
shows the deviation from the single-field prediction /nl = —5{ns — 1)/12 [27] represented by 
the green line in Fig. 9. Once again, multifield effects break this degeneracy. 
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Figure 8. Distributions for the amplitude of t lie power spectrum P^^, top left, scalar spectral index 
Us, top right, running of the spectral index, bottom left, and local non-gaussianity parameter /nl, 
bottom right. All distributions were taken for the sample with 2 active fields. 





Figure 9. Plot of the values of e at horizon-crossing versus r, left, and Ug versus /nl, I'ig'ht, for the 
sample with 2 active fields. The green lines represent the single-field predictions, r = IGe/Ta and 
/nl = -5{ns - 1)/12, respectively. 

5.3 Constraints from WMAP 

We now impose the observational constraints on the distributions. 

As can be seen in Fig. 10, the majority of the trajectories gives rise to values of running 
and spectral index that lie outside the observational bounds. From the running versus Ug plot 
alone, one could think that the trajectories giving rise to the peak around n^ = 1 would be 
in agreement with observations. Actually, a stronger constraint for n^ is imposed by the fact 
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Figure 10. Plot of the values of Ug versus r, left, and Ug versus running, right, for the model with 2 
active fields, blue points, and 6 active fields, red points. The lower panels greatly expand the vertical 
scales. The orange lines represent the 95% confidence limits using WMAP data. When running is 
allowed to be 7^ it is important to know what is the best scale with which to make comparisons 
with observations. In Ref. [26] this issue was explored; following their conclusions, we look at the 
constraints for k* = 0.017 Mpc^ which is not the same choice made in Ref. [25]. 

that the model predicts a negligible tensor-to-scalar ratio. Constraints on the n^ versus r plot 
shown in Fig. 10 exclude all the trajectories which result in n^ > 0.995. These constraints 
alone result in only ~ 10% of the total two-field sample, and ~ 50% of the six-field sample, 
of trajectories being in agreement with observations. 

A further constraint is imposed by requiring the correct amplitude of the scalar power 
spectrum, (2.5 it 0.1) x 10^^ [25]. Combining all constraints we obtained only two realisations 
in total concordance with observations in the full sample of 1086 cases of two-field inflation, 
and one realisation in the sample of 93 cases of six-field inflation. As discussed in the previous 
subsection, this is not a worrying result as the distribution of P^j; does not show a sharp peak. 

6 A close look at trajectories 

In this section we describe in detail the dynamics of individual trajectories and how they 
give rise to the distributions seen in the previous section. 

A remarkable feature of this model is that all inflationary trajectories encountered were 
essentially of the same type, inflection-point inflation with a wiggle in the trajectory. It 
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Figure 11. The dynamical betiaviour of Verse 20277, for tire region of the trajectory where inflation 
occurred. The top left plot shows the trajectory in field space plotted on top of the potential contours. 
The evolution of the amplitude of the power spectrum, top right, spectral index, bottom left, and 
running, bottom right, show different colours regarding the position in the potential; yellow represents 
evolution before the inflection point, orange between the inflection point and the fall of the ledge and 
red after the starting falling off the ledge. For the region represented in red, \ri\ > 1. 

turns out that while in principle inflation could occur at any location within the throat 
(above the tip), at least the last 60 e- folds of inflation always occur in a small sub-region 
(typically 0.02 < x < 0.09) in the vicinity of the inflection point discussed in §§2.2. While the 
inflection-point contribution to the potential, Vc + Vm, alone could not give enough inflation, 
the contributions from the bulk can alter the potential in such a way that sufficient inflation 
can occur. The range of contributions capable of giving inflation is limited so in practice the 
result is a very consistent dynamical behaviour. 

6.1 One inflationary trajectory to explain them all 

A typical inflationary realisation in our A'^p = 2 set-up is the affectionately-named Verse 
20277. Fig. 11 shows the inflationary trajectory superimposed on a contour plot of the 
potential, together with the evolution of the observable quantities we investigated. The 
trajectory evolves from right to left, passing an inflection point in the radial direction but 
also experiencing a slope in the angular direction, causing a wiggle in what would otherwise 
be standard single- field inflection-point inflation (see for example Ref. [28]). As will be shown, 
this behaviour explains the non-monotonic evolution of P(^(^, which we saw to be so common. 
Verse 20277, as ~ 84% of our realisations, has a blue spectral index, n^ > 1, which can also 
be understood by its inflationary dynamics. It also presents the most common evolution of 
/nl we encountered, giving rise to an unobservably small value of /nl, as seen consistently 
across our distributions. 
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Another interesting characteristic of the dynamics of Verse 20277, which again is repre- 
sentative of the whole sample, is that the process of falling off the ledge after the inflection 
point gives rise to a prolonged period of non-slow-roll inflation, often leading to 0(10) e-folds. 
This regime is shown in red in Fig. 11. While the slow-roll parameter e remains much smaller 
than 1 throughout all the analysis, the parameter \r]\ increases significantly, getting to values 
of ~ 35. The consequences of this effect will be discussed in what follows. 

6.2 A separable potential approximation 

A simple explanation for all of the above evolutionary traits presents itself if we approximate 
the inflationary region as a separable potential of the type.^ 



W{cP) = W{x, *) = C/(x) + 2^ Vii^i 

i=2 

where U{x) possesses an inflection point in the radial direction 



E 

i=2 



U{x) = Vi((/'i) = ao + ai{x - xq) + az{x 
and V{^) are slopes in the angular directions 



xq) 



(6.1) 



(6.2) 



(6.3) 



As already mentioned, for the majority of inflation, the evolution sits well within the slow-roll 
regime. During this period, we can write the number of e-folds of inflation as 



N{t^,U 



TVf 

E 



V,, 



(6.4) 



For this section we take the flat surface to be at horizon-crossing and the constant density 
surface to be at our time of evaluation. This allows us to write down the total derivative as 



dN 



^'E 



A^F 



E 



r*j 



'^j' 



(6.5) 



Following the procedure in Ref. [29] one can obtain an expression for the partial derivative 
in the previous expression as 




S: 



where the slow-roll parameters associated with each field are defined as 



1 



e,; = 



V. 



1 V 



m = 



,M 



i^^ 






2T3 \WJ ' ' T3W 
such that e = ^ Ej . Finally, we arrive at an expression for the derivatives of A'^ 



(6.6) 



(6.7) 



dN 



W* 



(6. 



Approximating small regions of potentials as separable is a very useful technique which can be applied 
quite generally. See Ref. [30] for another recent and rather nice example of this. 
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where we have the, rather important for the fohowing discussion, term 



^^^^ii'^^iA-^'h,) (6. 



which contains all the information about the constant density surface at the time of evalua- 
tion. All other terms are determined by horizon-crossing values. Note that due to assuming 
slow-roll, we have reduced our phase space to the N-p fields <f)i and that by approximating 
the potential as sum-separable, this expression no longer requires the 'c' and '*' surfaces to 
be infinitesimally separated. 

The second derivatives are given by [29] 

Having obtained expressions for the first and second derivatives of N, it is now possible 
to write down analytical expressions for the desired observable quantities. This is done by 
using the expressions from §§4.3 and using the horizon-crossing values for the correlation 
functions of 5(j). The amplitude of the power spectrum is given by 

nc-^^2^^ (6-11) 



where 



The spectral index is given by 



V* + Zj 



(6.12) 



1 _ Y^A^F Vt^f 
ns-l = -2e.-4 ^ J-' (6.13) 

Differentiating with respect to In k we also obtain a new expression for the running. We find 
it to be 



n' 



1 _ V ^^i^^ V 77*1/- (^ - ^^^^ V ''>■"? V ^> 
-8e*2 + 4Ve*r?*-16^ V^ - 8 ^ ^+4e* ^-2 

Z_-/ ' " / ..ax's ^-^ w? ^-^ uf 



1=1 



(6.14) 
Finally, the local non-Gaussianity is 



5 Z^j=l e* \^ 2e* j "•" Z^i,i=l e*e* ^J« 

^/nl = 2 ' ^ '^ ,,2 ^^ > (6-15) 

6 /y^Afp !^ 



^Here we have used the fact that d/dlnk « d/dN — {cj>i/H)d/d(j>i and continually made use of the 
substitutions N,iVi = -VT3 and <j>iN^ij/H = Vj/V + KijiY.i/VTg. 
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Figure 12. Example of evolution of the power spectrum Pc^q and the form of Zi giving rise to it for 
the separable potential. The power spectrum appears to stop evolving in the last few e-folds despite 
the bundle dilating in this region. 

where as a result of differentiating Zj, another term Aji containing contributions from the 
'c' surface is required [29]; 



dZ^ 




5r 



ik 



Vk' 



(6.16) 



All of the above expressions can be shown to reduce to the standard single-field formula by 
setting Ui = 1. 

With these expressions from the separable potential approximation, the phenomenology 
described at the beginning of this section becomes clear. We now discuss each of the trends 
encountered in the results individually. 

6.3 Why so many caterpillars? 

First, we would like to address the question of why non-monotonic evolution in the amplitude 
of the power spectrum is so common. 

Under the separable approximation, all superhorizon behaviour is encapsulated in the 
Zi terms. This must go to zero in the adiabatic limit since, as there is no evolution in this 
limit, the result must be independent of our choice of the time of evaluation. However this is 
not a sufficient condition for reaching the adiabatic limit; as we will see, it is possible for Zi 
to become negligible towards the end of inflation even when an adiabatic limit has not been 
reached. Before this limit is reached, Zi demonstrates significant variation. The extent to 
which Zi varies over the course of the last 55 e-folds of inflation changes between inflationary 
realisations according to both the realisation of the random coefficients and the path the 
inflationary trajectory takes in field space, but an example of the form of the variation is 
given in Fig. 12, where the fields have been redefined such that the inflection point is at the 
origin. 
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As the trajectory crosses the inflection point, Zj has a trough, then increases to zero as 
the trajectory leaves the plateau. As can be seen in Eq. (6.11), a peaking in \Zi\ corresponds 
to a peak in Pq(^, though the precise shape will also depend on the angular component which 
determines for how long the trajectory stays on the corresponding ridge in Zj. The peak in 
\Zi\ occurs due to e^ in the denominator of Eq. (6.9) reaching a minimum value and as such 
non-monotonic evolution is an inevitable consequence of the inflationary trajectory crossing 
an inflection point. The form of Zi also accounts for the final rise (the caterpillar's head) in 
the evolution of P(^^ since Zi turns up as a sum of quadratic terms in Eq. (6.11). This stage 
of the evolution is non-slow-roll but when comparing the evolution as given by Zi with our 
non-slow-roll transport code we found them to be in good agreement. 

6.4 Why so blue? 

Recognising that all inflationary trajectories take place in the vicinity of an inflection point 
also tells us about the spectral index. Remembering that Zj — )• towards the end of inflation, 
Eq. (6.13) states that the final value of the spectral index (along with the other observables we 
discuss) will be determined purely in terms of horizon-crossing values. This is what is often 
referred to as 'the horizon-crossing approximation' [31]. Unless very close to the inflection 
point, r]x dominates Eq. (6.13). When horizon-crossing takes place prior to the inflection 
point, rjx > 0, the horizon-crossing approximation says that n^ will necessarily be larger 
than 1. However if the situation arises where horizon-crossing takes place after the inflection 
point, r]x < and a red spectral index can be expected. This is what we find to be the case 
for approximately 16% of our trajectories (for Np = 2). The reason they are so rare is that 
60 e-folds of inflation still need to take place and this is difficult to achieve before falling off 
the ledge. It turns out that for the red cases, the total number of e-folds was always at least 
O(IOO). 

Fig. 5 shows a sensitivity to the number of fields A^f, where additional angular terms 
appear to act to redden the spectral index. This is not due to a larger proportion of trajecto- 
ries with horizon-crossing occurring after the inflection point. If this was the case then there 
would also be an increase in the number of e-folds. Instead, the increased redness seems to 
be a direct result of the geometry of the conifold and can be accounted for by the separable 
approximation. Writing Eq. (6.13) in a more suggestive form 

n,-l = -2ye* ^ + 2 '^ % , (6.17) 



i=l 






it is clear that additional angular terms act solely to redden the spectral index by suppressing 
the contribution from r]x as well as increasing e. 

6.5 Trends in the running 

Fig. 10 shows a lower bound in the plot of the running against the spectral index. This bound 
turns out to be parabolic and approximately proportional to — (n^ — 1)^, suggesting that there 
could be a dominant term in r;^ among the contributions for the running. Interestingly, in 
the single-field case, terms in rf do not contribute to the running, n' = — 24e^ + 16er/ — 2^. 
However, as can be seen by following Eq. (6.14), in the multifield case the situation changes. 
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Figure 13. Example of the form of Aji for the separable model. 



The terms from Eq. (6.14) in ry^. are: 
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which means that there is a residual r/^ cross term which dominates the expression. The 
parabohc lower bound is therefore intrinsically multifield. 

6.6 Evolution in /nl 

The separable approximation fares less well in reproducing the behaviour seen for /nl- The 
radial acceleration parameter ijx is the dominant contribution to Eq. (6.15), so the horizon- 
crossing approximation gives /nl < for horizon-crossing prior to the inflection point and 
/nl > for post inflection-point horizon-crossing. Bringing this together with the previous 
discussion for the spectral index implies there should be a correlation between Ug and /nl 
such that a blue spectral index is accompanied by /nl < and vice versa. This is in fact 
the relation from single-field inflation, /nl = —5{ns — 1)/12 [27]. While Fig. 9 clearly shows 
this line, it is also clear that the majority of points deviate from this trend. In fact, the 
histogram for /nl shows a tendency for /nl > 0, while the discussion until now would imply 
the opposite. 

Looking closely at particular evolutions of /nli like the one in Fig. 11, what we see is 
that during the non-slow-roll period at the end of inflation, the evolution of /nl experiences 
a rise, that typically forces the final value to be > 0. This behaviour is not predicted by the 
separable approximation. This discrepancy could be the result of modelling the inflection 
point as being cubic when in fact the ledge in the D-brane potential tends to be much more 
severe, or more interestingly it could be an intrinsically non-slow-roll effect. 

However, the slow-roll phase of the evolution does seem to be captured by the separable 
approximation. In addition to Zj, there is contribution to the evolution Aji which, as shown 
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Figure 14. Evolution of 0^^ top left, 9^ top right, and the dilation of the field bundle 9^^ bottom 
left for verse 20277. Bottom right is the histogram of the final field bundle widths for our two-field 
sample. The fields initially go through a region of focussing but all is undone as the trajectory falls 
off the ledge, resulting the persistence of isocurvature modes at the end of inflation. The histogram 
shows this to be true for all 1086 realisations, leading us to conclude the model is unpredictive without 
a description of reheating. 



in Fig. 13, has a peak and a trough. While the precise form varies between trajectories, this 
contribution gives rise to a peak and a trough in /nl, one above and one below the value 
/nl takes towards the end of inflation. 



7 How predictive? 

A consideration of paramount importance in any model of inflation is at what point does C 
cease to evolve. In single-field models, ( does not evolve on superhorizon scales but for any 
model with more than one field this tends not be the case. As we have already discussed, 
the criteria for whether or not evolution can persist is whether or not isocurvature modes are 
present. The limit in which these decay is referred to as the adiabatic limit. 

As already mentioned, one way to keep track of isocurvature modes is to monitor the 
dilation of the bundle of trajectories using Eq. (4.24). In the adiabatic limit the bundle is 
a caustic, for which a necessary but not sufficient condition is 9 — )• — oo, or equivalently 
— )• 0. Strictly speaking this limit cannot be reached during the slow-roll regime [17] but 
we can at least hope that isocurvature modes are exponentially suppressed. More generally 
^SR ^ Q corresponds to a region of focussing and 9^^ > to dilation, while 9^ > and 9^ < 
represent divergence and convergence to the momenta slow-roll attractor respectively. 
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It is interesting to note that our results show exceedingly consistent behaviour. As shown 
by the histogram in Fig. 14, in not one realisation did we find focussing to be occurring at the 
end of inflation. Rather, the process of falling off the ledge results in dilation of the bundle. 
This result is still present in the separable approximation and indeed we expect this result to 
apply to any model where inflation is terminated by falling off a ledge if sufficient focussing 
has not taken place previously. 

The downside of a lack of focussing at the end of inflation is that the model is technically 
unpredictive without knowledge of what takes place at the tip and the details of reheating. 
The plus side is that one should expect interesting evolution during subsequent periods. Since 
reheating is a non-linear process, this could give rise to interesting observational signatures. 
This is however beyond the scope of this paper. 

8 Summary 

In this paper we explored the multifield effects in D-brane inflation. To do this we made use 
of a particularly sophisticated model, originally developed in Ref. [2], to include contributions 
from the bulk containing random coefficients. These contributions make the multifield nature 
of D-brane inflation explicit. Our aim was to use the transport equations [3-5, 17] to study a 
large number of realisations of the potential, their resulting inflationary behaviours, and the 
consequent properties of C- We had two main objectives with this endeavour, first to present 
distributions for the observable predictions of this specific model and second to analyse how 
the isocurvature modes behave towards the end of inflation to understand how reliable those 
predictions are. 

Despite the random contributions, the inflationary behaviour was very consistent across 
different realisations. Inflation was always found to have canonical kinetic behaviour and to 
take place in the vicinity of the inflection point, which constitutes a very small region of field 
space. We also always found a break of the slow-roll condition |r/| < 1 when the infiationary 
trajectory falls off the ledge after the inflection point. This particular characteristic led to 
the extension of transport methods to non-slow-roll cases [5] . We also extended this method 
to compute the running of the spectral index. 

Regarding the predictions for observables, the amplitude of the power spectrum was 
found to be consistent with observation but unpredictive in so much as the spread in values 
was vast compared with the observational constraints. The spectral index tended to be blue 
but less so in the full six- field case. This fact, together with the negligible values found 
for the tensor-to-scalar ratio and the large spread in values of the spectral index running, 
moderately constrains the model given WMAP7 data. Without a reason to believe we are an 
atypical observer in our distribution, the soon-to-arrive data from Planck has the potential 
to put this model under considerable pressure. 

The non-Gaussianity parameter /nl was found to be consistently small. In over 1000 
inflationary realisations we found only 5 examples of potentially observable /nl- This fur- 
ther supports the recent discussions in Refs. [13, 30, 32] which suggests canonical multifield 
inflation is not expected to produce observable /nl- 

Due to the remarkable consistency in the inflationary behaviour, it was possible to ap- 
proximate the inflationary region of the potential as a sum-separable potential consisting of a 
polynomial inflection point in the radial direction and linear slopes in the angular directions. 
With this description all the observed characteristics of this model could be accounted for an- 
alytically except the final evolution in /nl • It was shown that the evolution of the amplitude 
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of the power spectrum undergoes non-monotonic evolution whenever the trajectory traverses 
the inflection point and this will typically give rise to a blue tilt which, along with all other 
observable quantities of interest, can be calculated purely in terms of horizon-crossing values. 
It was also shown that additional angular terms act to make the tilt more red. 

Regarding the reliability of these predictions, we made first steps towards a means of 
efficiently tracking isocurvature modes in non-slow-roll multifield models. Though we leave 
a complete description of this to a separate publication, we nevertheless found the bundle 
width to be exceedingly informative. Extending the ideas on the adiabatic limit discussed 
in Ref. [17] to non-slow roll constitutes a considerable increase in complexity. We found it 
helpful to consider two distinct bundles, one for the fields and a second for the momenta. 
With this we were able to conclude that in not a single realisation was an adiabatic limit 
reached. We expect this result to hold for any model where infiation is terminated by falling 
off a ledge, provided sufficient focussing does not occur prior to this. This result renders 
the model technically unpredictive, since subsequent events such as reheating can modify the 
characteristics of (^. 

Acknowledgments 

We would particularly like to thank David Seery for numerous discussions and guidance 
throughout this work. In addition we would also like to thank David Mulryne for advice 
on the transport method and Liam McAllister and Gang Xu for help with the building of 
the potential. M.D and J.F would also like to thank Tic Toe cafe, Brighton, where quite 
possibly the majority of this work was done. J.F also gives thanks to Leon Baruah for fixing 
his laptop at a critical stage of the writing up. M.D. was supported by FCT (Portugal), J.F. 
and A.R.L were supported by the Science and Technology Facilities Council [grant numbers 
ST/1506029/1, ST/F002858/1, and ST/1000976/1] and A.R.L. by a Royal Society-Wolfson 
Research Merit Award. 



References 

[1] D. H. Lyth and A. Riotto, Particle Physics Models of Inflation and the Cosmological Density 
Perturbation, Phys. Rep. 314 (1999) 1 [hep-ph/9807278]; D. Baumann and L. McAllister, 
Advances in Inflation in String Theory, Ann. Rev. Nucl. Part. Sci. 59 (2009) 67 
[arXiv:09Gl.G265]; D. H. Lyth and A. R. Liddle, The Primordial Density Perturbation, 
Cambridge University Press (2009); A. Mazumdar and J. Rocher, Particle physics models of 
inflation and curvaton scenarios, Phys. Rep. 497 (2011) 85 [arXiv: 1001.0993]. 

[2] N. Agarwal, R. Bean, L. McAllister and G. Xu, Universality in D-brane inflation, JCAP 09 
(2011) 002 [hep-th/1103.2775]. 

[3] D. Mulryne, D. Seery and D. Wesley, Moment Transport equations for non-Gaussianity, JCAP 
01 (2010) 024 [arXiv:0909.2256v3 [astro-ph.CO]]. 

[4] D. Mulryne, D. Seery and D. Wesley, Moment Transport equations for non-Gaussianity, JCAP 
04 (2011) 030 [arXiv:1008.3159v3 [astro-ph.CO]]. 

[5] M. Dias, J. Frazer, D. Mulryne, D. Seery and D. Wesley, Non-Slow-Roll Computation of C, by 
Transport: A How To, in preparation. 

[6] S. Kachru, R. Kallosh, A. D. Linde, J. Maldacena, L. McAllister and S. P. Trivedi, Towards 
Inflation m String Theory, JCAP 0310 (2003) 013 [arXiv:hep-th/0308055]. 



28 



\1 



C. p. Burgess, M. Majumdar, D. Noltc, F. Qucvcdo, G. Rajesh and R.-J. Zhang, The 
Inflationary Brane-Antibrane Universe, JHEP 0107 (2001) 047 [arXiv:hcp-th/0105204]. 



[8] D. Baumann, A. Dymarsky, S. Kachru, I. R. Klebanov and L. McAllister, D3-brane Potentials 
from Fluxes m AdS/CFT, JHEP 1006 (2010) 072 [arXiv:1001.5G28 [hep-th]]. 

[9] I. Klebanov and A. Murugan, Gauge/Gravity Duality and Warped Resolved Conifold, JHEP 03 
(2007) 042 [arXiv:liep-th/0701064]. 

[10] D. Baumann, A. Dymarsky, S. Kachru, I. R. Klebanov and L. McAllister, Holographic 
Systematics of D-brane Inflation, JHEP 0903 (2009) 093 [arXiv:0808.2811 [hep-th]]. 

[11] M. Tegmark, What does inflation really predict?, JCAP 0504 (2005) 001 
[arXiv:astro-ph/0410281v2]. 

[12] J. Frazer and A. R. Liddle, Exploring a String-Like Landscape JCAP 1102 (2011) 026 
[arXiv:1101.1619v2 [astro-ph.CO]]. 

[13] J. Frazer and A. R. Liddle, Multi-field inflation with random potentials: field dimension, 
feature scale and non-Gaussianity JCAP 02 (2012) 039 [arXiv:1111.6646vl [astro-ph.CO]]. 

[14] G. I. Rigopoulos and E. P. S. Shellard, Non-linear inflationary perturbations, JCAP 0510 
(2005) 006 [arXiv:astro-ph/0405185]. 

[15] G. I. Rigopoulos, E. P. S. Shellard and B. J. W. van Tent, Non-linear perturbations in 
multiple-field inflation, Phys.Rev. D 73 (2006) 083521 [arXiv:astro-ph/0504508]. 

[16] S. Yokoyania, T. Suyama and T. Tanaka, Primordial Non-Gaussianity in Multi-Scalar 
Slow-Roll Inflation, JCAP 0707 (2007) 013 [arXiv:0705.3178 [astro-ph]]; Primordial 
Non-Gaussianity in Multi-Scalar Inflation, Phys.Rev. D 77 (2008) 083511 [arXiv:0711.2920 
[astro-ph]]. 

[17] D. Scery, D. Mulryne, J. Frazer and R. Ribciro, Inflationary perturbation theory is geometrical 
optics in phase space, [arXiv:1203.2635vl [astro-ph.CO]]. 

[18] D. Secry and J. E. Lidsey, Primordial non-gaussianities from multiple-field inflation, JCAP 
0509 (2005) Oil [arXiv:astro-ph/0506056]. 

[19] G. Anderson, D. Mulryne and D. Seery, Transport method for the inflationary trispectrum, in 
preparation. 

[20] A. A. Starobinsky, Multicomponent de Sitter (Inflationary) Stages and the Generation of 
Perturbations, JETP Lett. 42 (1985) 152; D. H. Lyth, Large Scale Energy Density 
Perturbations and Inflation, Phys. Rev. D 31 (1985) 1792; M. Sasaki and E. D. Stewart, A 
general analytic formula for the spectral index of the density perturbations produced during 
inflation. Prog. Theor. Phys. 95 (1996) 71 [arXiv:astro-ph/9507001]; D. H. Lyth and Y. 
Rodrguez, The inflationary prediction for primordial non-gaussianity, Phys. Rev. Lett. 95 
(2005) 121302 [arXiv:astro-ph/0504045]. 

[21] M. Dias and D. Seery, Transport equations for the inflationary spectral index, Phys. Rev. D 85 
(2012) 043519 [arXiv:1111.6544vl [astro-ph.CO]]. 

[22] D. Wands, Local non-Gaussianity from inflation [arXiv:1004.0818vl [astro-ph.CO]]; D. 
Langois, F. Vernizzi and D. Wands Non-linear isocurvature perturbations and 
non-Gaussianities, JCAP 0812 (2008) 004 [arXiv:0809.4646v2 [astro-ph.CO]]; C. Byrnes and 
D. Wands Curvature and isocurvature perturbations from two-field inflation in a slow-roll 
expansion, Phys. Rev. D 74 (2006) 043529 [arXiv:astro-ph/0605679v3]. 

[23] D. Polarski and A. A. Starobinsky, Isocurvature Perturbations in Multiple Inflationary Models, 
Phys. Rev. D 50 (1994) 6123 [arXiv:astro-ph/9404061]; J. Garcia-Bclhdo and D. Wands Metric 
perturbations in two-field inflation, Phys.Rev. D 53 (1996) 5437 [arXiv:astro-ph/9511029]; D. 
Wands and J. Garcia-Bellido, Density Perturbations from Two-field Inflation, Helv. Phys. Acta 



29 



69 (1996) 211 [arXiv:astro-ph/9608042]. 

[24] E. Komatsu ct al., Seven- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 
Cosmological Interpretation, Astrophys. J. Suppl. 192 (2011) 18 [arXiv:1001.4538v3 
[astro-ph.CO]]. 

[25] E. Komatsu et al., Five- Year Wilkinson Microwave Anisotropy Probe (WMAP) Observations: 
Cosmological Interpretation, Astrophys. J. Suppl. 180 (2009) 330 [arXiv:0803.0547 
[astro-ph.CO]]. 

[26] M. Cortes, A. R. Liddlc and P. Mukherjee, On what scale .should inflationary observables be 
constrained?, Phys.Rev.D 75 (2007) 083520 [arXiv:astro-pli/0702170]. 

[27] J. Maldacena, Non-Gaussian features of primordial fluctuations in .single field inflationary 
models, JHEP 0305 (2003) 013 [arXiv:astro-ph/0210603]. 

[28] S. Hotchkiss, A. Mazumdar and S. Nadathur, Inflection point inflation: WMAP constraints and 
a solution to the fine-tuning problem, JCAP 1106 (2011) 002 [arXiv: 1101. 6046 [astro-ph.CO]]. 

[29] F. Vernizzi and D. Wands, Non-Gaussianities in two-field inflation, JCAP 0605 (2006) 019 
[arXiv:astro-ph/0603799v3], T. Battcfeld and R. Easther, Non-Gaussianities in Multi-field 
Inflation, JCAP 0703 (2007) 020 [arXiv:astro-ph/0610296]. 

[30] J. EUiston, D Mulrync, D. Sccry and R. Tavakol, Evolution of /nl io the adiabatic limit, 
[arXiv:1106.2153vl [astro-ph.CO]]. 

[31] S. A. Kim and A. R. Liddle, Primordial Non- Gaussianity in Multi-Scalar Inflation, Phys. Rev. 
D 74 (2006) 063522 [arXiv:astro-ph/0608186]. 

[32] J. Meyers and N. Sivanandam, Non-Gaussianities in Multield Inflation: Superhorizon 

Evolution, Adiabaticity, and the Fate of f^i^, Phys. Rev. D 83 (2011) 103517 [arXiv: 101 1.4934]. 



30 



